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Abstract 



^ ' OneLOop is a program to evaluate the one-loop scalar 1-point, 2-point, 3-point 

^ . and 4-point functions, for all kinematical configurations relevant for collider-physics, 
and for any non-positive imaginary parts of the internal squared masses. It deals 

■^ ■ with all UV and IR divergences within dimensional regularization. Furthermore, it 

[^ . provides routines to evaluate these functions using straightforward numerical inte- 

O ' gration. 
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1 Introduction 

The experiments at LHC for the study of elementary particles demand precise predictions of 
signals and potential backgrounds for new physics. This implies that for many processes the 
hard part has to be evaluated at least at the next-to-leading (NLO) level in perturbation theory, 
and this includes processes with four or more particles in the final state. There has been a 
considerable progress in this respect in recent years, resulting in full differential calculations 
for such processes [ni-fT2l. One of the bottlenecks that had to be resolved in order to achieve this 
was the calculation of the one-loop amplitude necessary to determine the virtual contribution, 
which in the mentioned publications involves at least 6-point functions. Developments have 
reached the stage that by now public programs exist that, given the machinery to calculate tree- 
level amplitudes, can calculate one-loop amplitudes involving 6-point functions and beyond at a 
high level of automation [|13[|14ll . 

Besides dealing with an increase in combinatorial complexity, moving from tree-level ampli- 
tudes to one-loop amplitudes demands the treatment of one-loop integrals. The most successful 
approach has been to express the one-loop amplitudes as a linear combination of universal one- 
loop functions, dependent on the particular scattering process under consideration only through 



the value of kinematical input parameters. The calculation of the one-loop functions is universal, 
while the calculation of their coefficients is not. In the "tensor-approach", these functions are the 
one-loop tensor functions. They involve a reduction scheme to express the tensor functions in 
terms of scalar functions, but this is universal. In the "unitarity-approach" the universal functions 
are the one-loop scalar functions themselves. 

As mentioned before, most approaches to calculate one-loop amplitudes demand the calcula- 
tion of the one-loop scalar functions. These may be infra-red (IR) divergent, and have to be dealt 
with within a regularization scheme. For NLO calculations it then suffices to have available the 
scalar 1-point, 2-point, 3-point and 4-point functions. Higher-point scalar functions can be ex- 
pressed as linear combination of these. The preferred regularization scheme for NLO calculations 
in QCD is dimensional regularization, and pubUc libraries LoopTools (VJ], QCDLoop [,1511 
and OneLOop [16J deliver all necessary scalar function for LHC kinematics and real particle 
masses. The first two rely on the public library FF for finite scalar functions IfTSl . 

Sooner or later, multi-particle NLO calculations will involve unstable particles. The preferred 
method to perform such calculations in a consistent way is the "complex-mass scheme", which 
demands the evaluation of one-loop functions with complex internal masses [[T9l . LoopTools 
provides those, including the most complicated finite 4-point function with 4 complex masses. 
The latter has been presented in [|20l , which is the completion and implementation of the work 
in [|2n . Alternative formulas for the 4-point function with complex masses have been given 
recently in [[22l|. 

In this write-up, the extension of OneLOop to deal with complex masses is presented. Es- 
sentially the only missing ingredient was the 4-point function with an arbitrary number of com- 
plex masses. It has now been implemented too, using the formulas from [|20l . As mentioned 
before, OneLOop provides all finite and IR-divergent scalar functions withing dimensional reg- 
ularization. It does not provide a complete treatment of exceptional phase space points corre- 
sponding to Landau singularities. Non- singular exceptional phase space points are taken care of 
to a certain extend, by expressing formulas in terms of differences of two dilogarithms divided by 
the difference of their arguments, and using numerically stable implementations of these. How- 
ever, as argued in [.22.1 . these exceptional cases do only rarely appear in practical calculations. 

Since the analytic continuation and the implementation of the formulas for the complex rou- 
tines is rather delicate, it is useful to have the means available to check their correctness nu- 
merically. First of all, they have successfully been compared with existing implementations. 
This includes a private code by A. Denner [|23l . As an extra check, they have been compared 
with results obtained via straightforward numerical integration. The methods applied to achieve 
this are presented in this write-up, and routines to evaluate all scalar functions involving at least 
two scales via numerical integration, thereby avoiding the evaluation of any dilogarithm, have 
been included in the OneLOop library. The actual integration is performed with the help of the 
CUBA-library [24| and routines from the Quadpack library [251 . 

OneLOop can be obtained from 



http : //helac-phegas .web . cern . chj 



2 Program description 

The main part of the package consists of 10 source files, written in FORTRAN 7 7 with common 
extensions, like enddo statements and do while statements. Furthermore, long names are 
used, and underscores are used in names. In particular, all names of subroutines, functions and 
common blocks begin with avh_olo_. It is in principle possible to compile the source files and 
link the object files directly to the main program the routines are supposed to be used in. It may, 
however, be more practical to create a library with the makefile coming with the package, and 
link the library instead. 

2.1 The 1 -point function 

The one-loop scalar 1 -point function is defined as 
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Ao = H(m., £) 



12 - m2 + ie ' 



(1) 



where the function h, of the renormalization scale \l and dimensional regularization parameter 
£ = 2 — D/2 follows the same convention as in ifTSl . and is given by 

P[l — £jr(l + £j ITT^ '^ 

The ie in the propagator denominators indicates the Feynman prescription for analytic continua- 
tion. Squared momenta follow the Bjorken-Drell convention l^ = Ig — if — l^ — if. 
The value of \x is by default equal to 1 . This default can be changed with 

subroutine avh_olo_mu_set ( mu ) 

The input mu shall be real, and has the dimension of energy, not energy squared. 
The 1 -point function is evaluated with 

subroutine avh_olo_aOc ( rslt ,m ) 

The input m shall be complex, and is the squared mass ra^ in Eq.O. It shall not have a positive 
imaginary part. If it does, an error message is returned, and the sign of the imaginary part is 
considered to be the opposite. 

The output rslt shall be a complex array of shape (0:2), with rslt ( ) containing the £° 
coefficient and rslt ( 1 ) containing the £^^ coefficient of the Laurent expansion in £ of Eq.©. 
The value of rslt (2 ) is alway equal zero. If m is equal zero, then rslt ( ) and rslt ( 1 ) 
are also zero. 



2.2 2-point functions 

The one-loop scalar 2-point function is defined as 
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where h,(|a., e) is defined in Eq.©. The 2-point function is evaluated with 

subroutine avh_olo_bOc ( rslt ,p,ml,m2 ) 

The input p, ml , m2 shall be complex, and is the square p^ of the momentum and the squared 
masses ra^ and ra2 in Eq.(|3]) respectively. The imaginary part of p shall be identically zero. 
If it is not, an error message is returned and the imaginary part is considered to be identically 
zero. The squared masses shall not have a positive imaginary part. If any of them does, an error 
message is returned, and the sign of the imaginary part is considered to be the opposite. 

The output rslt shall be a complex array of shape (0:2), with rslt ( ) containing the £° 
coefficient and rslt ( 1 ) containing the e^^ coefficient of the Laurent expansion in e of Eq.©. 
The value of r s It ( 2 ) is alway equal zero. If p and ml and m2 are equal zero, then rslt ( ) 
and r s It ( 1 ) are also zero. 

For the 2-point function, also the Passarino-Veltman functions Bi, Bqo and Bn are available. 
They are defined following 
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[F - m^ + ie] [(I + p)2 -mj + ie] 
They are evaluated with 

subroutine avh_olo_bllc ( bll, bOO, bl, bO ,p,ml,m2 ) 

The output bll,bOO,bl,bO respectively refers to Bi i , Boo and Bi and Bo. Each of them shall 
be a complex array of shape (0:2) again, with the entries containing the coefficients of the 
Laurent expansion in e. 



2.3 The 3-point function 

The one-loop scalar 3-point function is defined as 
P2 




P3= - Pi - P2 



Co = H(|j,, e) 
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where 



deni (I) = l^ — ra^ + ie 
den2(l) = [l + Vu 



vrq + ie 
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den3(l) = (l + pi +P2) -mj + ie , 
and H(|J., e) is defined in Eq.©. The 3-point function is evaluated with 
subroutine avh_olo_cOc ( rslt ,pl,p2,p3 ,ml,m2,m3 ) 
The input pi , p2 , p3 shall be complex, and in terms of the momenta in Eq.© it is given by 

Pi , P2 , vi = iV^+V2)'^ (8) 

The imaginary parts of these shall be identically zero. If any of them is not, an error message is 
returned, and it is considered to be identically zero. 

The input ml , m2 , mS shall also be complex, and in terms of the masses in Eq.© it is given 
by 

m^ , m^ , m^ (9) 

The imaginary parts of these shall not be positive. If any of them is, an error message is returned, 
and it will be considered to have the opposite sign. 

The output rslt shall be a complex array of shape (0:2), with rslt (0) containing 
the £° coefficient, rslt (1) containing the e^^ coefficient, and rslt (2) containing the e^^ 
coefficient of the Laurent expansion in e of Eq.®. 

By default, the values of r s It ( 1 ) and rslt ( 2 ) will be zero if the input does not exactly 
represent a configuration leading to an IR-divergent 3-point function. See Section [231 for more 
information. 



2.4 The 4-point function 

The one-loop scalar 4-point function is defined as 
P2 „ Ps 
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and H(|x, e) is defined in Eq.©. The 4-point function is evaluated with 



subroutine avh_olo_dOc ( rslt , pi, p2, p3, p4, pl2, p23 ,ml,m2,m3,m4 ) 

The input pi, p2,p3,p4,pl2,p23 shall be complex, and in terms of the momenta in Eq.(fTT]) 
it is given by 

Pi > P2 ) P3 > P4 = (Pl +P2 + P3)^ , (Pl+P2)^ , (P2+P3)^- (12) 

The imaginary parts of these shall be identically zero. If any of them is not, an error message is 
returned, and it is considered to be identically zero. 

The input ml , m2 , m3 , m4 shall also be complex, and in terms of the masses in Eq. (fTT]) it is 
given by 

mf , m2 , m3 , m4 . (13) 

The imaginary parts of these shall not be positive. If any of them is, an error message is returned, 
and it will be considered to have the opposite sign. 

The output rslt shall be a complex array of shape (0:2), with rslt (0) containing 
the e° coefficient, rslt (1) containing the e^^ coefficient, and rslt (2) containing the e^^ 
coefficient of the Laurent expansion in e of Eq.(fTOl). 

By default, the values of r s It ( 1 ) and rslt ( 2 ) will be zero if the input does not exactly 
represent a configuration leading to an IR-divergent 4-point function. See Section [23] for more 
information. 

2.5 IR-divergent one-loop functions 

OneLOop deals with all IR-divergent cases for the 3-point and the 4-point function within 
dimensional regularization. IR-finite cases are defined as those for which the e^ coefficient 
and the e^^ coefficient of the Laurent expansion in e are identically zero. IR-divergent cases 
are defined as those that are not IR-finite. They are all listed in [[T5l , and can also be found in 
Appendix 1X1 and Appendix |Bj 

By default, OneLOop only returns IR-divergent cases if the input is exactly referring to a 
IR-divergent case. This does not mean that the input has to be given in a certain order. Both 

zero = dcmplx ( OdO ) 

call avh_olo_cOc ( rslt ,m2,p2,m3 ,zero,m2,m3 ) 

and 

zero = dcmplx (OdO) 

call avh_olo_cOc ( rslt ,m3,m2,p2 ,m3,zero,m2 ) 

lead to a divergent case. The default does, however, mean that the input has to be numerically 
exact. For example, if an internal mass being equal to zero represents an IR-divergent case, then 
by default this divergent case is only returned if this mass is indeed identically zero on input. 
Also if an external momentum squared being equal to an internal squared mass represents an 
IR-divergent case, this divergent case is only returned if these numbers are indeed identical on 
input. The default can be changed with 



subroutine avh_olo_onshell ( thrs ) 

The input thrs shall be real, and is the threshold below which input is considered to be equal, 
or equal zero. It carries the unit of energy squared, so if for example the absolute difference of 
a squared momentum and a squared mass is smaller than thrs, they are considered to be equal. 
Also, if a squared mass or squared momentum is absolutely smaller than thrs, it is considered 
to be equal zero. If the user decides to stick to the default, warnings are returned whenever 
the input is such that it is close to a divergent case. The user is advised to take these warnings 
seriously. 

The scale associated with the regularization of the IR divergences is the same as the renor- 
malization scale |x of Section [ZTl 

2.6 Further useful routines 

All messages are by default sent to unit 6. This default can be changed with 

subroutine avh_olo_unit ( unit ) 

The input unit shall be integer, and if it is larger than 0, messages will be sent to unit unit. If 
it is not larger than 0, no messages will be printed at all. 

OneLOop can be ordered to print all input and output appearing during a run with 

subroutine avh_olo_printall ( unit ) 

The integer input unit is the unit to which all this information is sent. Again, if unit is not 
larger than 0, nothing will be printed. 

2.7 Real routines 

Some users may want to stick to real input. For all routines to evaluate one-loop functions, there 
are wrappers prepared that ask for strictly real input. For historical reasons, the routine names 
end on m rather than on r, and are given by 

subroutine avh_olo_aOm ( rslt ,m ) 

subroutine avh_olo_bOm ( rslt ,p,ml,m2 ) 

subroutine avh_olo_bllm ( bll,bOO,bl,bO ,p,ml,m2 ) 

subroutine avh_olo_cOm ( rslt ,pl,p2,p3 ,ml,m2,m3 ) 

subroutine avh_olo_dOm ( rslt , pi, p2, p3, p4, pl2, p23 ,ml,m2,m3,m4 ) 

2.8 Routines using numerical integration 

As mentioned in the introduction, with the OneLOop package also come routines to evaluate 
one-loop functions using numerical integration. These have been included as means of a cross- 
check. The main purpose is to check the correctness of analytic continuations, and the emphasis 
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has been put on straightforwardness rather than sophistication in their implementation. In par- 
ticular the routines for the 4-point function can be slow, and will give only the first few decimals 
correctly. Since errors in analytic continuation lead to rather drastic differences, not just in the 
last few decimals, this is in general sufficient. 

The routines come as an independent library, which from user point of view essentially is a 
copy of the main OneLOop library. It has all routines, except those for the 2-point Passarino- 
Veltman functions, with the only difference that the routine names start with avh_oni instead 
of avh_olo. It has to be compiled separately, and needs access to the cure routine of the 
Cuba library [f24l . It provides two extra routines, namely 

subroutine avh_oni_maxeval_set ( maxeval ) 
subroutine avh_oni_unitcuba_set ( unit ) 

The integer input maxeval of the first routine changes the maximal number of integrand eval- 
uations cure is ordered to perform, which by default is 10''. The integer input unit of the 
second routine changes the unit to which messages related to the numerical integration are sent, 
which is 6 by default. If this number is not larger than 0, no messages will be printed. 

2.9 Quadruple precision 

Some compilers provide the option to compile at a higher precision level, such that for example 
all 16-byte complex variables become 32-byte variables. All series expansions inside OneLOop 
have been prepared to also deal with this precision level. During its first execution, OneLOop 
will determine the precision level at which it has been compiled, and from then on use the relevant 
expansions. 

3 One-loop scalar functions via numerical integration 

In this section, integral representations for the one-loop scalar functions suitable for numerical 
integration are derived. Since the aim is to have an alternative numerical implementation of 
the scalar functions in order to check the correctness of the analytic continuations in the "hard- 
wired" implementation, the main guideline in the derivation is that it should be straightforward. 
Methods to smoothen the integrand in order to speed up numerical integration, as for example 
used in [I 26ll27l , are not applied. Other approaches to numerical integration of loop integrals can 
be found in [l28H30l[3Ml. 

The main problem in the numerical integration of loop integrals is the existence of singulari- 
ties in the integrand. These can be categorized into two groups: those that can be dealt with via 
deformation of the integration contour, and those that cannot. The former will be dealt with. . . via 
contour deformation. Among the latter are the so-called mass-singularities, or infra-red (IR) sin- 
gularities. They cause the integral to be divergent, and are typically inherent to the particular 
scattering process under consideration. They ask for a regularization scheme, and are dealt with 



within dimensional regularization in this section. Other singularities that cannot be dealt with 
via contour deformation typically appear in isolated points in phase space. A complete treatment 
of those, like for example in [27|, is beyond the scope of this work. 

In this section, we write the one-loop n-point scalar integral within dimensional regulariza- 
tion in the following convention 
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(14) 



The denominator momenta q, are related to the external momenta pj via c\^ = and 

j-i 

k=1 



(15) 



The reason for the choice of the factor in front of the integral is that it leads to a simple expression 
for the scalar function in terms of the so-called Feynman parameters, namely 
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where S is the linear operator represented by the symmetric u>i n matrix with elements 

I] L; -r I] Li, — I u; — uv 

Sj,k 



(16) 



m? + m^ 



1) 



qkj 



The integration region is such that all variables Xj run from to cxd. 

On = [0, oof, 



(17) 



(18) 



however the Dirac delta-function in Eq.(fT6l) confines the integration region eventually to the n- 
simplex. One could proceed by direct elimination of one of the integration variables with the 
help of the delta-function. We shall take a slightly different strategy, and scale the variables xi 
to Xn-i by Xn, which we call ij, to find 
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Here, we denote 
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We used the fact that we are allowed to replace e by y^e since ij^ is positive. Elimination of the 
y -integral now gives 

(l+&i^i)^""" 
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(21) 



The integration region is not finite anymore. One advantage of expressing the one-loop integral as 
Eq.(l2T]) is that one does not have to bother about the end-point at infinity if one wants to explicitly 
deform the integration contour in the complex plane in order to avoid poles as prescribed by ie. 
It does not matter towards which direction in the complex plane infinity is approached, since the 
integrand does not have a singularity at infinity. Another advantage is that a proliferation of terms 
is prevented when one wants to apply techniques involving integration by parts, since integrands 
will be at infinity. 

3.1 Contour deformation 

The ie -prescription dictates that the contour deformation has to be such that the imaginary part 
of the polynomial in the denominator of the integrand in Eq. (|2T1) is negative when the real part 
vanishes. Furthermore, the deformation has to be such that the end-points at Xj =0 (the begin- 
points) stay fixed. IR-singularities, which cannot be dealt with by contour deformation, may 
show up as vanishing real parts at Xj = 0. 

Since we want to allow for squared masses with non-zero, but negative, imaginary parts, the 
matrix A and the vector b may have imaginary parts. We re-define our notation, and write 

Af-A-ir , b^b-iy, (22) 

now with A and b real. The matrix V and the vector y have only non-negative components. 
Our deformation is very similar to the one in [32]. We introduce a strictly increasing function 
T such that t(0) = and t(oo) is finite, for example t(x) = TooX/(1 + x) with Too > 0. The 
deformation is given by 

X t-^ z{x] — X — iy[x] y (23) 

with 

yj(x)=T(Xj)(Ax + b)j, (24) 

and where (Ax + b)j denotes the j-th component of the vector Ax + b. The Jacobian matrix 
coming with the deformation is given by 

9xk 
The imaginary part of the polynomial in the denominator of the integrand in Eq. (l2T)) satisfies 



1 -iT'(xj)(Ax + b)j 5j,k-iT(xj)Aj,k. (25) 



Im 



|z-(A-ir)z + 2(b-iy)-z+cj = 2y-(Ax + b) -y-Ty +x-rx + 2y-x-Im{c} . (26) 



The last three terms are safe, in the sense that they never become negative. The first term never 
becomes negative if y is chosen as in Eq. (|24)) (notice that the function t acts as a deformation of 
the inner product, keeping it however positive definite). We have to worry a bit about the second 
term. The matrix V has elements 

n.k = Tlj + rik , Tij = -llm{m- } > , (27) 
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so 

'n-1 \ /n-1 



Applying the Schwartz inequality to both sums, we find 



n-l 



|y-ry|<2v/^^^^||Tl||||g||^ , ||Tif = }^Tif. (29) 



Now we have for the first two terms on the r.h.s. of Eq. (|26l) 

^-^^^ 

and we conclude that we are safe of we choose t such that 

-1/2 
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2y-iAx + b]-y-ry>2'^i^-V^^^^\mvf] , (30) 
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[n-1)}^Im{mf}M . (31) 

j=i / 



Many methods for numerical integration ask for integrands living on the unit hypercube. For 
our application, this can conveniently be achieved by preceding the contour deformation with the 
mapping 

[0, 1] ^ [0, oo) , pj ^ Xj = 1/pj - 1 , (32) 

leading to a Jacobian weight factor pr^ = (1 + Xj]^ for each integration variable Pj in the 
integrand. 

3.2 IR-divergent scalar integrals 

For IR-finite integrals, the integrand in Eq.(l2T1) can be expanded in e and integrated term by 
term. In particular, e can be put to to calculate the £°-term. For IR-divergent integrals, this 
is not directly possible. In this work, we are only interested in scalar integrals up to n = 4, 
i.e. box integrals and triangle integrals. IR-divergent bubble integrals vanish within dimensional 
regularization. 

Regarding the box integrals, we use the fact that the IR-divergent structure of one-loop in- 
tegrals with more than 3 external lines can be expressed in terms of triangle integrals. This 
is completely worked out in [|35l . Formulated for our application, that work explains how to 
determine the coefficients Cj for every IR-divergent box integral such that it can be written as 

4 

U = ^c^hii)-lT (33) 
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where I3 ( j ) is the triangle integral obtained from the box integral by removing propagator de- 
nominator j, and where I4" is IR- finite 111 Putting the terms under one integral and with a common 
denominator, we can write 



If = ^Cjl30]-l4 
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(35) 



j=i j=i j=i 

The integral in Eq.(|34l) is IR-finite, and can be integrated numerically. The appearance of the 
term y in the numerator seems counter-intuitive in this respect, and in fact the coefficients Cj 
appear to be such that it vanishes for every IR-divergent box integral. In terms of Feynman 
parameters, I4" is given by 



rfin 

-■■4 



d^xSl 1 

04 






(x(x-Qx— ^i^x-Sx) — p-x + y 



fx-Sx — ie 



|2+£ 



(36) 



where the matrix Q and the vector (3 have components 



(37) 



To summarize, the IR-divergent box integrals are calculated following Eq. (|33l) . where I4" is 
finite and can be integrated numerically directly. The coefficients Cj for the triangle integrals 
can be determined as explained in [35J, and are explicitly given in ^!I\. They can also be found 
in Appendix lAl together with the coefficients a and |3j. For the IR-divergent triangle integrals 
themselves integral representations are derived in Appendix |B] which can be performed numeri- 
cally. 



4 Summary 

The program OneLOop for the evaluation of one-loop scalar functions up to and including 4- 
point functions was presented. It deals with all kinematical configurations relevant for collider- 
physics, and any non-positive imaginary parts for the internal squared masses, and it deals with 
all UV and IR divergences within dimensional regularization. Furthermore, it provides routines 
to evaluate all these functions using straightforward numerical integration, and the methods ap- 
plied to achieve this were presented. The program can be obtained from 



|http : /7helac-phegas .web . cern . ch/ 



' Within the convention of Eq.(fT4b. I3 ( j ) actually is a triangle integral divided by — ( 1 + £ ) 
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A Coefficients for IR-divergent box integrals 



In this appendix, the coefficients Cj for Eq. (|33l) and a, (3j for Eq.(l36l) are given. The coefficient 
y in the latter equation vanishes for all cases. The coefficients are expressed in terms of the 
external momenta defined as in Eq.(fTOl). and we denote 



Sl2 = (Pl +P2)' 



S23 = (P2+P3]' 



(38) 



For completeness, we also give the matrix S in Eq. (|36l 



/ 



1 

2 






V 



mt 



2mf 
fm| 
hm? 



mi 



Pt 

S12 

■vi 



raf + ra2 - 
2m| 
vnl + vnl- 
ral + m| - 



Pf 



Vl 



mf 
ml 



$23 mf 



-ml- 

fmf ■ 

2mi 

fm^ ■ 



S12 

vl 



vl 



mf 

m^ 

mf 



fm|- 

-m| - 

fm|- 

2m| 



^ \ 



P4 

S23 

vi 



and the matrix Q in Eq.(l36l) 
/O 



1 

2 



v 





2pf 

Sl2+Pf - 
Pf+P|- 







\ 



S12 



p^ 



Pi 



HPt-P2 Pl- 
2si2 S12 

^'-vi 



-vl- 
-vl- 



S23 
-vl 



J 



(39) 



(40) 



S23 S12 + P4 

The numbering of the IR-divergent box integrals refers to the numbering in [[T5l . Some of the 
separately labeled box integrals in that publication are just particular cases of other box integrals 
with one or more propagator masses equal to zero. Only the general cases are drawn below. The 
graphs follow the same conventions as in fTSl]. Dashed lines correspond to vanishing invariants. 
Explicitly mentioned invariants are non-zero and not equal to other explicitly mentioned invari- 
ants. External lines with the same color as internal lines, and with the same label in the form 
of a squared mass, indicate that the squared external momentum is equal to the squared internal 
mass. 



Box (16) and (15): one soft and no coUinear singularity. ra| may vanish. 



Pi 



mi 



Pl 



m. 



m 



/ \ 

mi (16) ml 



1 



S12 - mf 
1 



a 



(0,0,1,0) 



S12 - mf 

P = 2 (^ ' ^12 - P2 + ^2 ) 2si2 , S12 - S3 + w}^ 

S12 - mf ' 



(41) 
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Box (14): two soft and no coUinear singularities. 



m| 



\ 



m 



/ 



/ 



m 



\ 



m: 



mi (14) m^ 



Sl2 



:i, 0,1,0] a=— P = (0, 1,2,1). 

Sl2 



(42) 



Box (13), (12), (11), (10), (9) and (5): one coUinear singularity, and no, one, or two soft 
singularities, ra^ and ra4 may vanish in box (13), and ra4 may vanish in box (12). 



pi 2 P3 



m' 



\. 



m; 



\. 



P^ 



m' 



m 



m 



\ , 

(13) pI 



(12) 



"\ 



Pi 




A = (si2-m3)(s23-m4] - (p2-m3)(p4-ra4 
c = ^(0, 0, S23-P4, S12-P2) 



oc — 



S12 + S23 -Vi-vl 



1 



2 2 
4-P3 



|3 = - (0 , , si2(si2 + 2s23 -Vi-Vi- Vl) + P2(P3 " vl) 



, P4(Si2 + S23 + P3 - P4 - ^Vl) + S23(Sl2 " P3)] 



(43) 



Box (8), (7), (6) and (4): two coUinear and at least one soft singularity. ra4 may vanish in 
box (8). 



./ 



pI 



/ 



m 



./ 



m 



m 



(8) 



■\ 



Pi 



m 



(7) 



\ 



P^ 



m 



\ 



(6) m 
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1 



a 



|3 



Sl2lS23-Ta^J 

S12 + 2S23-P4-P3 



(S23-P3> 0, S23-P4, S12) 



Sl2lS23-TTL|j 



1 



sn[S23-Tnl 



[0,0, suisn + 2s23 - P4 - P3) 

, P4(Sl2 + S23 + P3 - Vi) + S23(Sl2 " P3)) • 



(44) 



Box (3), (2) and (1): more than two coUinear, or two coUinear and no soft singularities. 

\ '' ^ '' 



Pi 



(3) 



■\ 



P'4 



(2) 



■\ 



P^4 



(1) 



c = 



S12S23-P2P4 



/ 2 2 2 2\ 
(,S23 — P2 > S12 — P4 , S23 — P4 , S12 — P2J 



a = 2 



S12 + S23-P2-P4 

2 2 

S12S23 ~ViVa 



P = a(0, 0, S12, P4) . 



(45) 



B IR-divergent triangle integrals 

In this appendix the integral representations for the IR-divergent triangles are given, suitable 
for numerical integration. The aim is to arrive at those with as few algebraic manipulations 
as possible, starting from the representation of Eq.(l2T]). This will be done only for triangles 
involving at least two mass scales (besides the extra scale \i due to the regularization), since 
the analytic formulas for the triangles involving only one mass scale are rather simple. The 
considered triangles are depicted in FigUl following the labeling as in ifTSl . The general triangle 
modified Caley matrix is 



(46) 



1 j 

2 1 


f 2mf 


2 1 2 


-vi 


mf + mf - 


-vl 


mf + ml - 


- Pf 2mi 




ml + vnl- 


-vl 


V mf + mf - 


-pf m^ + m|- 


-Vi 


2ml 
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m 




Figure 1: IR-divergent scalar triangles with at least two mass scales, labeled as in ifTSl . Triangle 
(2) in that publication is a particular case of triangle (3) with ra^ = 0. 



Triangle (6) 

Triangle (6) has modified Caley matrix 



1 








2m| 



ra| + ra| — P2 



2mi 



The integral representation Eq. (|2TI) of this triangle is given by 



roo 



13(6) 



roo 



dx 



^-l+2e 



dij 



(m|y^ + (m| + m| - p|)ij + mj 



2ll+e 



The x-integral can be performed directly, leading to 



13(6) = / 

2£ 



dy 



;i+y 



|2£ 



(m^y^ + [mj + ml- p^y + miy+' 



(47) 



(48) 



(49) 



The remaining integral is convergent for any value of e. The integrand can be expanded into a 
Taylor series in £, which can be integrated numerically term by term. 



Triangle (3) 

Triangle (3) has modified Caley matrix 



1 

2 

















m^- 


vl 
















m^- 


vl 


m' 


— 


vi 


m' 


— 


Pi 


2Tn 


1 



(50) 



We choose to eliminate row and column 1 instead of 3 to obtain the integral representation 
Eq. ((2T]) . leading to 



l3(3) 



roo 



roo 



dx 



dy 



;i +x+-y; 



-1+2e 



mV + (m^ -pf]xy + (m^ -p|)x]i+« ' 



(51) 
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Whereas the integral for triangle (6) was divergent due to the large-x behavior, this integral is di- 
vergent due to the small-x behavior. Notice that a factor x^ +^ can be taken out of the denominator. 
Using this fact, we apply subtraction to write 



l3(3) 



dx 



dy 



;i+x + y] 



-l+2e 



X 



l+£ 



1 



1 



((m2-p2)y + m2-p2) 



2^1- 



+ 



n^/HTi-e) 



poo 



w 



dy 



:i 



\-l+e 



((m2-p2]y + m2-p|] 



2^^+^ ' 



where we used the identity 



roc (1+^+ )-i+2. 7rVM-r(-£) ^^ , i+, 
o'^ ^^ ^ r(l-e) ^^+^^ 



(52) 



(53) 



Both integrals in Eq.(l52]) are now finite. To obtain the £^-term and the £°-term, the integrand 
of the second integral may be expanded to order £, and e may be put to zero in the first integral. 
Notice that the first, 2-dimensional, integral vanishes for ra^ = 0, i.e. for triangle (2) in IfTSl . 



Triangle (4) 

The integral representation for triangle (4) becomes 



l3(4) 



poo 



poo 



dx 



dij 



:i +x + y 



1-1 +2£ 



Rescaling the variable y <— xij, we get 



l3(4) 



poo poo 



dx 



dy 



and performing the x-integral, we find 

l3(4)==l 



.m^x^ + [vn^ - Tpl)xyy+^ ' 
(1 +x + xy]-^+^' 

(1+y)'^ 



Substituting y <— 1 /ij we get 

Is (4) 



-1 



dy 



(1+y)'^ 



1 



,l+£ 



[m^y + m^ - Vly+' ' 



(54) 



(55) 



(56) 



(57) 



which can be dealt with using subtraction 

l3(4) 



2£ 



Mr. (1+y)^- 
dy r 

y^+^ 



1 



1 



+ 



7ti/M^rf-£i 



m^y + m^-pjy+'' 
1 



:i+y)(m2-p^ 



2 _T,2ll+£ 



r(i-£) [m^-viV^' 



(58) 



The remaining integral is finite, and can be expanded in £ on the integrand level. 
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